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0\ . Abstract 

^ ' An Ansatz for the curvature contribution to the density of states for massive 
1—5 ! 

OO I quarks in a bag is given and shown to reproduce exact mode-fihing calcula- 
tions. A mass- formula for spherical lumps of 3- flavor quark matter is derived 

>. 

,—1- I self-consistently from an asymptotic expansion within the MIT bag model, 

CO ' taking into account bulk, surface, and curvature contributions. Shell model 



calculations are performed for a variety of strange quark masses and bag con- 



o 

■^ ' stants, and the results shown to match nicely with the asymptotic expansion. 
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Several experiments with relativistic heavy- ion coUiders [Q, as well as cosmic ray searches 
have been performed or planned to test the possible (meta) stability of strange quark matter; 
a phase of 3-flavor quark matter, that could be the ground state of hadronic matter (thus 
more stable than iron) even at zero temperature and pressure. The first notion of this 
possibility seems to be due to Bodmer [Q], and much work, also on the consequences for 
neutron stars and cosmology, has been performed since the idea was resurrected by Witten 
1^. (See Ref. [Q for reviews and references on strange quark matter). 

It is therefore important to know the signatures to be expected for small lumps of quark 
matter (strangelets) composed of up, down and strange quarks. An important property 
is the mass (or energy per baryon). A calculation involving mode-filling in a spherical 
MIT-bag was performed for ud-sjsteras by Vasak, Greiner and Neise 0, and for 2- 
and 3-flavor systems for a few parameter sets by Farhi and Jaffe [0 and Greiner et al. [§. 
Recently, Gilson and Jaffe |^ published a thorough investigation of low-mass strangelets for 
4 different combinations of s-quark mass and bag constant. Such calculations, capable of 
showing shell-effects etc., are rather tedious, but important for an understanding of decay 
modes. For many applications, including generalization to finite temperature, a global mass- 
formula analogous to the liquid drop model for nuclei is of great use. 

An investigation of the strangelet mass-formula within the MIT bag model was performed 
by Berger and Jaffe 0]. That investigation included Coulomb corrections and surface 



tension effects stemming from the depletion in the surface density of states due to the mass 
of the strange quark. Both effects were treated as perturbations added to a bulk solution 
with the surface contribution derived from a multiple reflection expansion. 

Recently it was pointed out that another contribution to the energy, the curvature term. 



is dominant (and strongly destabilizing) at baryon numbers below a hundred |[TT|. A mass- 
formula (liquid drop model) for systems of massless quarks including the curvature was 
derived in [|T2| and shown to fit well with shell model calculations. 

However, the density of states correction due to curvature was only known for mass- 
less quarks (and for infinite mass quarks, where the problem corresponds to the Dirichlet 



boundary conditions studied by Balian and Bloch |T3|), whereas the surface tension is a 
function of the mass, vanishing for zero quark mass. Thus low-mass u and d quarks could 
be consistently included, but for realistic s-quark masses (expected to be in the range of 
100-300MeV) neither m^ = nor uis — ^ cxd are terribly good approximations. It was argued 
in []rT| that it seemed as if the shell model calculations of Ref. [^ required roughly 2 mass- 
less curvature contributions plus the surface energy. As will be shown below, the s-quark 
curvature contribution indeed goes through zero for intermediate s-quark mass. 

In the ideal Fermi-gas approximation the energy of a system composed of quark flavors 
i is given by 

E = Y,{n, + N,^jii)+BV (1) 

i 

Here VLi, Ni and //« denote thermodynamic potentials, total number of quarks, and chemical 
potentials, respectively. B is the bag constant, V is the bag volume. 



In the multiple reflection expansion framework of Balian and Bloch [0, the thermody- 
namical quantities can be derived from a density of states of the form 

""'■ A'S.^um^s.fomcX, (2, 
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where area S = § dS {= AtxR^ for a sphere) and curvature C = § i^ + -^ j dS (= 87ri? for 
a sphere). Curvature radii are denoted Ri and i?2- For a spherical system Ri = R2 = R. 
The functions fs and fc will be discussed below. 

In terms of volume-, surface-, and curvature-densities, riiy, rii^s, and rii^c, the number 
of quarks of flavor i is 

f'<^Ft dNi 

Ni= —rr^^ = ^i,vV + ni,sS + Ui^cC, (3) 

JO dk 

with Fermi momentum kpi = {jJ^I — mfY^'^ = fii{l — XfY^"^; \i = rrii/ jjii. 
The corresponding thermodynamical potentials are related by 

a = a,yV^ + a,5^ + a,cC, (4) 



where dfli/dfii = —Ni, and dflij/dfii = —Uij. The volume terms are given by 
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The surface contribution from massive quarks is derived from fs{m/k) 
{1 - (2/7r)tan-i(A;/m)}/87r as [|10| 
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For massless quarks Vti^s = '^i.s = 0, whereas /c(0) = —1/2471^ gives ||TTl , p!2| Vti^c = /^i/S^r^; 
«i,c = -fii/Ai:^- 

The curvature terms have never been derived for massive quarks, but as will be shown 
in this Letter, the following Ansatz (found from analogies with the surface term and other 
known cases) works: 

1 
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This expression has the right limit for massless quarks (/c = —1/2477^) and for infinite 
mass, which corresponds to the Dirichlet boundary conditions studied by Balian and Bloch 
ITSl (/c = 1/127?^). Furthermore, the expression gives perfect fits to mode-filling calcula- 
tions (see the Figures and discussion below). From this Ansatz one derives the following 
thermodynamical potential and density: 
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With these prescriptions the differential of E(y, S, C, Ni) is given by 

dE = J2 i^iydV + Vti^sdS + Vti^cdC + ^ndNi) + BdV. (12) 

i 

Minimizing the total energy at fixed N^ by taking dE = for a sphere gives the pressure 
equilibrium constraint 

B = -Y. ^.V - I E ^^,S - J^ E ^^,C. (13) 

i i i 

Eliminating B from Eq. (^ then gives the energy for a spherical quark lump as 

E = Y.{N,fi, + \q,,sS + ^a,cC). (14) 

The importance of the curvature contribution is demonstrated in Figure |l| which shows 
the energy per baryon for fixed bag constant of a bag filled with one flavor of quarks, quark 
mass ranging from 1 to 450 MeV. The shell model results are derived in a similar manner 
as described by Gilson and Jaffe [^], except that no zero-point correction has been added. 
One sees that a significant positive curvature energy arises for low quark mass (where the 
surface term approaches zero). For a quark mass around 150 MeV (for the given choice 
of B) the curvature term is zero, and the surface term alone reproduces the shell model 
calculations, whereas for higher rrii the surface term overshoots, but the shell model results 
are reproduced when including the now negative curvature contribution. For high masses one 
approaches the Dirichlet boundary conditions, but since quark masses are always lower than 
the chemical potential rather than infinite, the Dirichlet-limit is only a rough approximation, 
and the full surface and curvature terms should be used. 

Equally nice fits are obtained for other choices of B and mg. Small changes in the 
Ansatz for the curvature density of states leads to clearly visible errors in the fits to the 
mode-filling calculations, so even though the curvature term has not been derived from 
the multiple reflection expansion (and attempts to do it have been unsuccessful), it seems 
fair to conclude, that this is the proper term, and that surface plus curvature corrections 
are sufficient (and equally necessary!) to reproduce the overall behavior of mode-filling 
calculations. 



Figures ^ and |^ compare shell model and liquid drop model calculations for three-flavor 
strangelets for massless u and d quarks, and massive s quarks for a variety of parameters. 
Again, the overall behavior of the energy per baryon is fitted perfectly by the liquid drop 
results. For low rris, the energy is roughly that of bulk plus 3 massless curvature contribu- 
tions. For intermediate mg, 2 massless curvature terms plus the s-quark surface term are 
important, and for high s-quark mass, where no s-quarks are present, the energy per baryon 
is given by E/A in bulk plus the curvature terms from 2 massless quarks. (The approach of 
Ref. |jlO| gives constant E/A in the low as well as high ms-limit since only surface tension 



is included, and for intermediate values of rus the significant contribution of u and d quark 
curvature is lacking). 

The calculations neglected Coulomb energies, the phenomenological zero-point energy, 
and the gluon exchange described by the strong coupling, a^. Coulomb energies can be 
self-consistently included (see |]12|); but are at most a few MeV. The zero-point energy was 
included in as — 1.84/i?, and reduces the energies for very small A somewhat. In E/A 
this term contributes (for fixed chemical potentials, which is not quite correct for A — » 1 
| |14]| ) like A~^^^, compared to curvature (A^^/^) and surface {A~^^^). Significant changes 
occur only for A < 10, and can also be accounted for by an extension of the liquid drop 
model. There is no easy way to include the strong coupling. 

A detailed study of strangelet decay modes (an important issue for experimental 
strangelet searches) must ultimately rely on shell model calculations. The main scope of this 
Letter has been to show, that the general structure of shell model results for multi-quark 
systems can be explained by a liquid drop model mass-formula including volume, surface, 
and curvature contributions. Since such calculations are much easier to perform than the 
shell model calculations, this allows the overall structure of strangelets to be studied for a 
range of parameters, and facilitates generalizations to finite temperature, studies of phase 
transitions and mixed phase structures in heavy-ion physics, cosmology and neutron stars 
T3[ . For a consistent liquid drop treatment, the full curvature term for a massive quark. 



introduced here for the first time, turns out to be crucial. 
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FIGURES 
FIG. 1. Energy per baryon of 1-flavor quark-lumps at fixed bag constant {B^''^ = 145MeV), 
for quark masses of 1, 150, 300 and 450 MeV (bottom to top). For each mass 4 curves are shown: 
Horizontal lines are the bulk values; dotted curves bulk plus surface terms; spiky curves are the shell 
model calculations, and smooth curves on top of the shell model results are the liquid drop model 
results including surface and curvature. Notice that the dotted curve is almost indistinguishable 
from the bulk curve for mg = IMeV, and overlaps with the full calculation for nis = ISOMeV. 

FIG. 2. Energy per baryon for 3 flavor strangelets for B^' = 145MeV (the lowest possible 
value for ud-quark matter to remain unstable) . Pairs of curves with the same signature show shell 
model and liquid drop model calculations for the same nis ranging from 50 to 300 MeV in steps of 
50MeV (bottom to top). For even higher nis no s-quarks are present, and results look like those 
for nis = 300MeV. 

FIG. 3. As Figure 2, but for B^''^ = 165MeV (the highest possible value for stable uds-quaik 
matter with nis = 0). Now the highest nis that can be distinguished is 350 MeV. 
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